Solar system constraints on the Dvali-Gabadadze-Porrati braneworld theory of gravity 
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A number of proposals have been put forward to account for the observed accelerating expansion 
of the Universe through modifications of gravity. One specific scenario, Dvali-Gabadadze-Porrati 
(DGP) gravity, gives rise to a potentially observable anomaly in the solar system: all planets would 
exhibit a common anomalous precession, dui/dt, in excess of the prediction of General Relativity. 
We have used the Planetary Ephemeris Program (PEP) along with planetary radar and radio 
tracking data to set a constraint of \dcu/dt\ obs < 0.02 arcseconds per century on the presence of any 
such common precession. This sensitivity falls short of that needed to detect the estimated universal 



precession of \dui/dt\ DGP = 5 x 10~ 4 arcseconds per century expected in the DGP scenario. We 
discuss the fact that ranging data between objects that orbit in a common plane cannot constrain 
the DGP scenario. It is only through the relative inclinations of the planetary orbital planes that 
solar system ranging data have sensitivity to the DGP-like effect of universal precession. In addition, 
we illustrate the importance of performing a numerical evaluation of the sensitivity of the data set 
and model to any perturbative precession. 
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fcH ! I. INTRODUCTION 

OJj. 

Measurements of the apparent brightnesses of Type la supernovae provide strong evidence that the Universe is 
accelerating in its expansion [l], Q. Although General Relativity (GR) with a non-zero cosmological constant can 
accommodate these observations, it is a major challenge for theoretical physics to explain the observed value of 
the energy density of the vacuum [H, 0] • This cosmological constant problem has motivated the theoretical physics 
community to explore alternative descriptions of gravity that can account for the accelerating Universe without the 
need for dark energy. 

ly-^ , One class of such theories includes the five-dimensional braneworld gravity developed by Dvali, Gabadadze and 
Porrati (DGP gravity) in which standard model interactions are constrained to a four-dimensional brane, embedded 
in a higher-dimensional bulk. The graviton is also constrained to the brane on small scales, but over cosmological 
■ distances is free to explore the full higher-dimensional space, thereby weakening gravity on Gigaparsec length scales. 
The appealing feature of this theory is its ability to produce cosmic acceleration without invoking dark energy. The 
one free parameter in DGP gravity is the crossover scale r c , above which gravity becomes five-dimensional. In order 
to produce the observed acceleration rate of the Universe, r c rss 5 Gpc. 

The self-accelerating phase of DGP gravity has been shown to contain a ghost at the linearized level (fj. This 
suggests that the self- accelerated phase may be unstable and therefore not self-consistent. For the purposes of this 
work, we point out that the interpretation of the ghost mode is a subject of debate in the current literature @, [8], and 
that there is at least one way to preserve self-acceleration while avoiding the ghost instability Q. We also note that 
the "normal" (non-accelerating) phase of DGP gravity, which does not possess a ghost @, also predicts a periapse 
precession of the same magnitude but opposite sign from the self-accelerating phase and is thus constrained by 
the current work. 

Lue and Starkman [l(| derived the remarkable result that on the scale of the solar system, DGP gravity predicts a 
universal rate of periapse precession for bodies in nearly circular orbits, independent of not only the orbital scale but 
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also the mass M of the central body: 

3c / r c \~ l 

= =5x 10~ 4 ( - — ) arcseconds per century (1) 

8r c \5 Gpc / 

where c is the speed of light. Thus, upper limits on the rates of precession of the planetary perihelia set a lower 
limit on r c , the scale above which the DGP graviton is free to explore the bulk. To see that the precession rate is 
independent of the orbital semi-major axis and of M, consider the weak-held potential of any gravitational theory 
V(r) = Vjv(r) + SV(r), where Vjv(r) = —GM/r is the Newtonian potential, and 6V(r) is a perturbation such that e = 

5V(r)/Vfj{r) « 1. For nearly circular orbits, the precession per orbit Slu can be computed from Su> = ixr [r 2 (r _1 e)'J 
[ilj . where ' = d/dr. The precession rate dw/dt is found by dividing 5w by the orbital period P. For a Keplerian 
orbit, P = 27rr 3 /VGA/. In DGP gravity, if r << r* then SV(r) = cVGMr/r c [loj], which leads to a precession rate 

1/3 

that is independent of both r and M. We have made use of a new scale, = vr g rt) , below which the metric is 
dominated by the central body rather than by the effects of the cosmological flow [10l ], Here, r g = 2GM/c 2 is the 
gravitational radius of a mass M. For the Sun, M = 2 x 10 30 kg, r g 3 km, and so r* ss 130 pc. Therefore all solar 
system bodies easily satisfy r << r*. Iorio [Hj] extended Equation [j] to second order in the orbital eccentricity e. He 
showed that e does not appear at first order and so for solar system bodies, the eccentricity contributes less than 1% 
to dw/dt. We ignore the contribution of the eccentricity to the precession in this work. 

Using solar system ephemeris constraints from p]| , Iorio has undertaken a program [l2|, EH, EH EH to constrain 
anomalous gravitational effects using solar system observations. Our approach here is similar, but with three dif- 
ferences: (1) We use an independent ephemeris code, the Planetary Ephemeris Program (PEP) [l7| to arrive 
at limits on anomalous precession. (2) We undertake a constrained global fit to the motions in the solar system, 
as advocated by Nordvedt [l8j]. (3) We perform a simulation to quantitatively assess our sensitivity to a DGP-like 
precession. To test the DGP scenario, we constrain the anomalous precession rates of the planets (beyond the GR 
prediction) to be equal, as specified by Equation [1] 

The recent renaissance of Lunar Laser Ranging (LLR) , brought about by the millimeter- pre cision range measure- 
ments from the Apache Point Observatory Lunar Laser-ranging Operation (APOLLO, see [lj|), promises improved 
sensitivity to lunar orbital parameters. As discussed in [111 ]. LLR observations could provide substantially improved 
tests of the DGP scenario. As seen in Equation [1] DGP theory predicts the same amount of periapse motion for the 
Moon as for the planets, even though they orbit different central bodies. We will defer a consideration of improved 
lunar ranging limits on DGP to a future publication. Our focus here is to explore DGP limits based on ranging 
between planets in orbit about the Sun. 



~dt 



II. THE PEP EPHEMERIS MODEL: 
USING PRECESSIONS TO LINK OBSERVATIONS WITH THEORY 

PEP integrates the Parametrized Post-Newtonian (PPN) equations of motion for massive bodies [2(| HH, with 
free parameters that describe the masses and orbital properties of solar system objects. The equations of motion 
also encode deviations from standard physics that include, but are not limited to, the violation of the equivalence 
principle and the time rate of change of the gravitational constant. The model parameter estimates are refined by 
minimizing the residual differences, in a weighted least-squares sense, between observations and model predictions. 
This procedure can be iterated until convergence criteria are met. 

The PEP software has enabled constraints on departures from standard physics. For example, it has been used to 
place limits on the PPN parameters (3 and 7 (see e.g. [22|, HH) and on possible violations of Lorentz Invariance [24| . 
For the results described here we fixed the PPN parameters to their GR values (j3 = 7 = 1), and we modified the 
code to allow for a common perihelion precession for all planets in excess of the GR predictions. This approach allows 
us to quantitatively assess the covariance between the anomalous precessions and the more mundane solar system 
parameters such as planet masses and initial orbital elements. 

Our modification to the PEP code allows for independent limits on any anomalous precession of the planets in the 
solar system, or more generally, for constrained ratios of precession rates for the different solar system objects. Rather 
than generating a specific version of PEP to encode the equations of motion for each speculative gravity scenario, we 
favor taking this perturbative approach based on anomalous precessions. For the results described here we constrained 
all precession rates to be equal, as predicted by DGP gravity. 

For a broad class of extensions to GR, the PPN formalism 0, [2l[ has been a very useful framework. Unfortunately 
many of the interesting ideas that have been put forward recently (motivated by the Dark Energy crisis) do not readily 
map onto the PPN parametrization. We are not aware of an analog of the PPN formulation that readily encompasses 
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the diverse modified-gravity scenarios that are currently under consideration. For example, even a simple Yukawa 
modification to GR does not naturally lend itself to a PPN description. The "forward modelling" approach, to look for 
the (necessarily small) corrections to GR dynamics in the solar system in the form of perihelion precessions, provides 
a convenient tie point between observation and a broad class of speculative theories. 

In fact, a bound, two-body system will exhibit closed orbits for only two central potentials: V(r) oc r 2 (e.g. Hooke's 
law), and V(r) oc 1/r (e.g. Newtonian gravity) pEI ]. A perturbation to the Newtonian 1/r potential will, in general, 
violate Kepler's First Law by generating unclosed orbits. These orbits are characterized by the rate of precession of 
the argument of periapse (or perihelion when the body orbits the Sun). The argument of periapse u>, defined in the 
orbital plane, is the angle between the ascending node and periapse (closest approach). The argument of periapse 
should not be confused with the longitude of periapse, w, which is the sum of u> and the longitude of the ascending 
node f2. In general, the angles u> and fl are non-coplanar and so w is often called the broken angle. In Figure [TJ 
L02 = ^2 — CI, and w\ — w\ (the longitude of the ascending node for the orbit of mass 1 is undefined). In a purely 
Kcplerian orbit, ui, CI and w do not change over time. 

Perihelion precession has been used before to distinguish theories of gravity. For example, GR garnered early success 
from its ability to explain the long-standing 43 arcsecond per century discrepancy between the observed advance of 
Mercury's perihelion and the Newtonian prediction. This difference is caused by the 1/r 3 perturbation to the effective 
gravitational potential in the weak- field, slow- motion limit of GR [26| . Similarly, any modification to gravity or 
new fundamental interaction that deviates from a 1/r potential could be tested through its effect on the periapse of 
solar system bodies. A useful framework for computing the expected precessions for a broad class of non-l/r central 
potentials has been developed by Adkins and McDonnell [271 ] ■ We invite proponents of modified-gravity models to 
provide the expected ratios of periapse precessions in the solar system, and we would be pleased to undertake a PEP 
analysis at these fixed precession ratios. 



III. DGP PRECESSION, GLOBAL ROTATION AND THE NEED FOR INCLINED ORBITS 

One might legitimately worry that a set of differential observations taken within the solar system, namely the 
ranges between the Earth and other solar system bodies, might not be sensitive to a universal precession of the kind 
predicted by DGP gravity. As we discuss below, differences in the inclinations of the planetary orbits with respect 
to the ecliptic are vital for testing DGP in the limit of small eccentricity, as they break the degeneracy between a 
DGP-like precession and a global rotation. 

As described by Williams and Standish (28|, range data between two objects in heliocentric orbits can be used to 
determine the two orbital semi-major axes and the difference between the orbital angular frequencies of the bodies. 
If the orbits are coplanar, then the ranging data are not sensitive to a common (DGP-like) precession of the perihelia 
of the two orbits, because the constant excess rotation rate vanishes from the orbital angular frequency difference. In 
other words, a rotation preserves both the lengths and the relative angle between vectors. Adding additional coplanar 
orbits and associated range data does not help. Absent a determination of the phase of one of the objects in its orbit 
around the Sun, a DGP-like common precession would remain undetectable. For non-coplanar orbits, however, range 
data between planets can distinguish a global rotation from a DGP-like precession, given adequate measurement 
accuracy. 

To see that non-coplanar orbits are required to distinguish DGP gravity from a global rotation, consider the range, 
pit), between the centers of mass of two planets in orbit about the Sun with position vectors r±(t) and r^if) 

P 2 (t)=r 2 1 +r 2 2 -2r 1 -r 2 (2) 

where r = |r| and the labels 1 and 2 refer to the reference body (e.g., the Earth) and the body on the inclined orbit 
(the ranged-to body), respectively. In an elliptical Keplerian orbit, r is related to the angle from periapse (called the 
true anomaly), /, by 



1 + e cos / 

where e is the eccentricity of the orbit, and a is the orbital semi-major axis. From the law of cosines in spherical 
trigonometry, the angle between the position vectors of the two bodies is 

= cos(A2 — CI) cos(Ai — CI) + sin(A9 — Cl) sin(Ai — Cl) cosi. 

nr 2 

where Ai = w\ + /1 and A2 = vd 2 + f% are the true longitudes of the bodies, vo\ and -072 are the longitudes of periapse 
for the two bodies, Cl is the longitude of the ascending node and i is the mutual orbital inclination. Figure [1] depicts 
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the relevant geometry. In an unperturbed Keplerian orbit, w, and i are static, and only /, the angular position of 
the body, measured with respect to periapse, changes with time. Periapse precession causes w to change in time as 
well. 

Following Williams and Standish [28[, we expand Equation^ to first order in the orbital eccentricities and to second 
order in the inclination (the inclination appears only in even powers) 

p 2 (t) = a\ + a 2 2 ~ 2a x a 2 cos(A 2 — Ai) + E + I (3) 

where 



E = — 2a 2 ei cos/i — 2a 2 e 2 cos / 2 
+ 2a\a,2 [ei cos /i cos(A 2 — Ai) 
+ e 2 cos/ 2 cos(A 2 - Ai)] 



and 

/ = S i n 2 i [ C0S (A 2 - AO - cos(A 1 + A 2 - 20)] . 

To first order in the eccentricity and second order in the inclination, Equation [3] is equivalent to Equation 3 of 28;] , 
which is written in terms of mean anomalies and longitudes. The first three terms in Equation [3] simply give the 
range between circular, coplanar orbits. The function E accounts for the effect of the eccentricities of the orbits, and 
I arises from the relative inclination of the orbits. It is immediately clear that the longitudes appear only through 
the difference, A 2 — Ai, except in /. Therefore measurements of p(t) cannot distinguish between a DGP-like universal 
precession and a net rotation unless i 7^ 0. Furthermore, the sin 2 i term strongly suppresses the effect of a DGP-like 
precession on the observable range to planets whose orbits are only mildly inclined to the orbital plane of the reference 
planet (e.g., the ecliptic). 



IV. DGP PERTURBATION TO PLANETARY RANGES 



Given that the planetary orbits (save for the Earth's) in our solar system are inclined to the ecliptic, we can estimate 
the size of the DGP perturbation to interplanetary ranges. Over a span of time, At, the total DGP-induced perihelion 
precession is Aui = (dui/dt)At << 1. To first order in Aw, the difference, Ap(t), between the DGP range prediction 
and the Newtonian range prediction, pN(t) given in Equation is 

Ap(t) = -^ 2 -Awsin 2 zsin(A 1 + A 2 - 20). (4) 
2p N (t) 

In our solar system, the planetary orbital inclinations with respect to the ecliptic are small, and therefore the sin 2 i 
factor strongly suppresses the DGP signature. For example, the orbits of Mercury and Mars are inclined to the ecliptic 
by 7.0 and 1.8 degrees, and so sin 2 i — 10~ 2 and 1CP 3 , respectively. Unfortunately, there are no range measurements 
of Pluto, whose orbital inclination is 17 degrees (and so sin 2 i « 0.1). Using Equation[4j we have computed Ap(t) over 
40 years for Mercury and Mars using the DGP prediction of \dio/dt\ — 5 x 10~ 4 arcseconds per century. The range 
differences, shown in Figure [H do not exceed one meter. As described in Section IVII the single-point measurement 
uncertainties in the ranging data to Mercury and Mars are approximately 200 m and 5-40 m, respectively. 



V. OTHER SOURCES OF PRECESSION 



The planetary perihelia precess from Newtonian perturbations alone, by amounts much larger than relativistic 
effects. For example, Newtonian perturbations from other planets induce a perihelion precession in Mercury's orbit 
of more than 500 arcseconds per century. In addition, the Newtonian potential of the oblate Sun deviates from that 
of a spherically symmetric body, causing the perihelia of the planets to precess. A planetary orbit that is inclined by 
an angle 6 to the rotation axis of the Sun will experience the following rate of perihelion precession due to the Sun's 
mass quadrupole moment [26| : 

dLj ) =%( *° y^r-^w^- (5) 



dtJj 2@ P Wl-e 2 ) 



FIG. 1: Two planets, Mi and M2 (assumed massless), in Keplerian orbits about a common focus (e.g., the Sun). The perihelia 
are labeled pi and p2- The longitudes of perihelion, vji and zz>2 are measured with respect to the reference direction x, with 
W2 the algebraic sum of non-coplanar angles. The true anomalies, /1 and /2, are measured with respect to their respective 
perihelia, and the true longitudes are defined by Ai = zu\ + fi and A2 = T&2 + ji- The inclined orbit has two additional 
parameters, the longitude of the ascending node, Q, and the inclination, i. In our solar system, the planetary orbits are inclined 
by only a few degrees with respect to the ecliptic (the plane of the Earth's orbit about the Sun). 

Here a is the orbital semi-major axis, e is the orbital eccentricity, P is the orbital period of the planet, Rq is the mean 
radius of the Sun, and J20 is the dimensionless coefficient of the quadrupole moment of the Sun's mass distribution. 
For comparison, the classic GR prediction for the perihelion precession of an orbiting planet is 



where Mq is the mass of the Sun. This term, which has a different radial dependence than does the effect of the solar 
oblateness, produces the predicted 43 arcseconds per century perihelion motion for Mercury mentioned earlier. The 
amounts of perihelion precession for each planet contributed by Equations [5] and [5] are listed in Table fl] 



There is a rich, publicly available archive of range and Doppler measurements to the terrestrial planets, spacecraft 
and landers. Table HT1 summarizes the data that were used in this analysis, including the type and duration of the 
observations and the typical associated measurement uncertainties. 

The range and Doppler measurements to Mercury and Venus were obtained by reflecting radar signals from the 
planets, and therefore depend on the planetary surface topography. Similar radar measurements of Mars exist but 




(6) 



VI. OBSERVATIONS 
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(a)Mercury (b)Mars 

FIG. 2: The estimated evolution of the range difference, Ap(t), between DGP and Newtonian gravity to the planets Mercury 
(left) and Mars (right). The difference of zero at the start of the time span reflects the freedom of the model to choose the 
initial anomaly / of each planet independently. 



TABLE I: The perihelion precession rates, in arcseconds per century, from GR and the solar oblateness, computed from 
Equations [5] and [6] with Mq = 2 x 10 30 kg, Rq — 6.95 x 10 5 km, J20 — 2x 10 -7 , and solar equatorial orbits (8 = tv/2) assumed 
for all bodies. The DGP prediction for r c — 5 Gpc is 5xl0~ 4 arcseconds per century, which is less than the GR effect for 
r < 37 AU, but exceeds the precession from the oblate Sun for planets outside the Earth's orbit. As in the text, a is the orbital 
semi-major axis (AU), and P is the orbital period (years). 



Planet 


a 


P 


( — ) 


1 dui \ 

V dt I GR 


Mercury 


0.39 


0.24 


2.5 x IQ-" 


43.0 


Venus 


0.72 


0.62 


2.6 x 10~ 3 


8.6 


Earth 


1.00 


1.00 


8.4 x 10" 4 


3.8 


Mars 


1.52 


1.88 


1.9 x 10~ 4 


1.4 


Jupiter 


5.20 


11.9 


2.6 x 10" 6 


6.2 x 10" 2 


Saturn 


9.54 


29.5 


3.1 x 10~ 7 


1.4 x 10~ 2 


Uranus 


19.2 


84.0 


2.7 x 10 -8 


2.4 x 10 -3 


Neptune 


30.1 


164.8 


5.6 x 10" 9 


7.8 x 10~ 4 


Pluto 


39.5 


248.5 


2.4 x 10 -9 


4.2 x 10~ 4 



have been supplanted by superior radio tracking to Mars orbiters and landers, most notably the Mars Global Surveyor 
(MGS) and Mars Odyssey orbiters. Range observations to the Jovian planets are extremely limited. There are four 
radio tracking observations of spacecraft as they passed by Jupiter (Pioneer 10 and 11 and Voyager 1 and 2) and two 
such observations for Saturn (Voyager 1 and 2). No observations of Uranus, Neptune or Pluto were used here. More 
detailed information about the planetary range and Doppler data sets is provided by Standish [29(. LLR data were 
not used in this work. This allows for a more direct comparison with the previous results of Pitjeva [13j |. and the 
subsequent DGP interpretations by Iorio [3, EH, EH • 

VII. RESULTS 

Using the input data described above, and allowing for a single common anomalous precession for the planets 
Mercury, Venus, Earth, Mars, Jupiter and Saturn, the PEP code derived a universal precession rate of (dui/dt) bs = 
0.001 ± 0.002 arcseconds per century (la statistical), consistent with no anomalous precession. In Section IVIII1 
however, we show that this statistical uncertainty overestimates the sensitivity to precession by a factor of 11. Over 35 
years (the approximate duration of our data set) a precession rate of 0.002 arcseconds per century would have perturbed 
the orbit of Mercury and Mars by «100 m and ~500 m, respectively, leading to observable range perturbations, Ap(t), 
of ~2 m and «1 m. For comparison, an independent analysis of planetary data by Pitjeva [13j, in which each planet 
was allowed to precess independently, derived la statistical uncertainties in the estimates of the planetary perihelia 
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TABLE II: Characteristics of the data used in this analysis. N t s and a obs are the number of observations and typical single- 
point measurement uncertainty for each data set. For range measurements, the uncertainty is given in meters of one-way path, 
while for Doppler data the fractional frequency uncertainty Ai//v is reported. MGS stands for the Mars Global Surveyor. 
Fly-by refers to spacecraft that passed close to Jupiter and Saturn (see text for details). For reference, the orbital semi-major 
axes of the planets are given in Table [I] The line-of-sight velocities of Mercury and Venus, as seen from the Earth, can reach 
35 km s _1 and 15 km s^ 1 , respectively. 



Planet 


Type 


Target 


Date Range 


N obs 


17 


obs (m or ^) 


Mercury 


Radar range 


Surface 


07/1969 - 08/1997 


7,880 




200 




Radar Doppler 


Surface 


07/1969 - 03/1974 


174 


(5 


- 20) x 10~ n 


Venus 


Radar range 


Surface 


07/1969 - 07/1982 


5,466 




200 




Radar Doppler 


Surface 


07/1969 - 05/1977 


340 


(0.2 


- 20) x 10 -11 


Mars 


Radio tracking 


Mariner 9 


11/1971 - 10/1972 


185 




500 




Radio tracking 


Viking 1 & 2 


07/1976 - 11/1982 


3,872 




15-40 




Radio tracking 


Pathfinder 


07/1997 - 09/1997 


90 




40 




Radio tracking 


MGS 


02/1999 - 06/2005 


124,856 




12 




Radio tracking 


Odyssey 


02/2002 - 06/2005 


142,416 




6 


Jupiter 


Radio tracking 


Fly-by 


1973, 1974, 1979 


4 




4,500 


Saturn 


Radio tracking 


Fly-by 


1980, 1981 


2 




4,500 



of 5.0, 300, 0.4 and 0.5 milliarcseconds per century for Mercury, Venus, Earth and Mars, respectively. Thus Pitjeva's 
tightest constraint on planetary precession, allowing for the independent motion of each body, is five times smaller 
than our statistical constraint on the uniform precession rate. We raise a flag of caution regarding the use of published 
solar system precession constraints based on statistical errors only. In particular, we show in the next section that 
the statistical uncertainty of 0.002 arcseconds per century overestimates the true sensitivity of PEP and the data set 
to a DGP-like precession by a factor of 11. 



VIII. ESTIMATING THE SENSITIVITY TO DGP PRECESSION 



In order to determine the sensitivity of PEP and the data set to a universal, DGP-like perihelion precession, a 
series of simulation data sets were generated with the characteristics of our real data (see Table [TTJ) , but that also 
had varying amounts of excess perihelion precession, ranging from zero to 0.2 arcseconds per century, the latter being 
100 times the statistical limit reported in the previous section. The resulting best-fit DGP-like universal precession 
rates, as determined by PEP, were consistently 11 times smaller than the input precession rates. This exercise showed 
that ~90% of the precession signal was absorbed into other ( "nuisance" ) model parameters due to the non-linear 
dependence of the model on the anomalous precession rate. To account for this dilution, we choose to report an 
uncertainty on (cko/dt) b s that is equal to the amount of input precession required to change the central value of that 
parameter by one statistical standard deviation (in this case 0.002 arcseconds per century, as reported in the previous 
section). Our realistic la constraint on a common DGP-like precession is therefore: 



dw 
~dT 



< 0.02 arcseconds per century. 



IX. IMPLICATIONS FOR DGP GRAVITY 



The universal precession rate for the DGP scenario depends only on r c , the scale above which the graviton can 
escape the four-dimensional brane and explore the full five-dimensional space. In order for DGP gravity to explain 
the observed acceleration rate of the Universe, r c must be ~5 Gpc. Our ler upper limit on the common anomalous 
planetary precessions in the solar system of \dw/dt\ < 0.02 arcseconds per century requires that r c > 0.13 Gpc. This 
is not yet stringent enough to rule out the DGP scenario as being responsible for the accelerating expansion. 



X. CONCLUSION 



Constraints on the periapse precession of solar system bodies can be used to compare solar system observations 
with speculative gravitational theories, especially those that fall outside of the PPN framework. We have modified 
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PEP to allow for anomalous precessions of the planetary perihelia, and have performed a least-squares solution of 
this modified model to a suite of archival radar and radio tracking data, with the requirement that all anomalous 
precession rates be equal. At the level of 0.02 arcseconds per century, we find no evidence for a universal precession 
in excess of the GR prediction. This result, which places a lower limit of 0.13 Gpc on the cross-over scale in DGP 
gravity, is a factor of ~ 40 away from the prediction of DGP gravity. We therefore cannot rule out DGP gravity as 
a mechanism for the accelerating Universe. It is remarkable, however, that a theory of gravity designed to explain 
phenomena on cosmological scales has the potential to be tested within the solar system. 

We have pointed out that mutually inclined orbits are necessary to distinguish a universal DGP-like precession from 
a global rotation, and that the observable signature of DGP gravity for interplanetary ranges scales like sin 2 i. We 
also showed that the statistical uncertainty on the precession rate, determined by a weighted least-squares fit to the 
observations, overestimates the sensitivity to a DGP-like precession by a factor of 11. One should therefore exercise 
caution when interpreting statistical uncertainties from the literature as constraints on physical theories. 

This work can be extended to incorporate the nearly 40 years of high-precision lunar laser ranging observations, 
including the new, millimeter-precision, LLR data from the APOLLO project. In addition, the recent successful asyn- 
chronous laser transponder [30[ experiment to the MESSENGER spacecraft over 24 million km [3l[ has demonstrated 
that ~20 cm precision planetary ranges over interplanetary distances are possible. Currently, a feasibility study is 
underway for a project to deploy a laser transponder on the surface of Mars capable of producing millimeter-precision 
range measurements [32j . Such a high- precision planetary range data set, if extended over several years, would greatly 
improve the solar system constraints on gravitational physics, and could potentially rule out the DGP scenario. 

To facilitate the tie between observation and theory, we encourage others to calculate the relative precession rates 
of solar system bodies in speculative gravitational theories that do not have a PPN expansion. 
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